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Role of GTP remnants in microtubule dynamics 
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We study a model of microtubule assembly /disassembly in which GTP bound to tubulins within 
the microtubule undergoes stochastic hydrolysis. In contrast to models that only consider a cap of 
GTP-bound tubulin, stochastic hydrolysis allows OTP-bound tubulin remnants to exist within the 
microtubule. We find that these buried GTP remnants enable an alternative rescue mechanism, and 
enhances fluctuations of filament lengths. Our results also show that in the presence of remnants 
microtubule dynamics can be regulated by changing the depolymerisation rate. 
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Microtubules are semiflcxible polymers that serve as 
structural components inside the eukaryotic cell and are 
involved in many cellular processes such as mitosis, cy- 
tokinesis and vesicular transport [H, H, In order to 
perform these functions, microtubules (MTs) continually 
rearrange through a process known as dynamic instability 
(DI), in which they switch from a phase of slow elonga- 
tion to rapid shortening (catastrophe), and from rapid 
shortening to growth (rescue) [ij. Despite extensive ex- 
perimental efforts over several decades, it has not been 
possible to discriminate between models that have been 
proposed to explain DI d, H, 0, Q. In this letter, 
we study a minimal model of DI that involves stochas- 
tic hydrolysis (SH), a mechanism that has received rela- 
tively little attention compared to interfacial hydrolysis 
(IH) that forms the basis of cap models H. With re- 
cent advances in experimental techniques0, fiot . it has 
become possible to quantify MT dynamics at nano-scale 
and, thereby, provide more stringent tests of models. We 
make predictions that clearly distinguish between SH and 
IH mechanisms, and experimental tests of these predic- 
tions should clarify the nature of hydrolysis in MT dy- 
namics. 

MTs are formed by assembly of a — /3 tubulin dimcrs, 
which are polar and impart polarity to MTs. MTs grow 
mainly from the end that has exposed /3 tubulin, and are 
composed of (typically) 13 linear protofilaments.[l| 

While a free tubulin dimer has a GTP molecule bound 
to each monomer, incorporation into a MT activates the 
/3-tubulin monomer for hydrolysis of its associated GTP. 
GDP-bound tubulin is less stable within the MT lattice 
[m and hence a GDP-bound tubulin at the tip of a MT 
has a higher rate of detachment (depolymerization) than 
a GTP-bound tubulin. While it is established that GTP 
hydrolysis is essential to DI, the mechanism of the process 
is not fully understood. Models in the IH class assume 
that all hydrolysis occurs at a sharp interface between 
GDP-bound and GTP-boimd tubulins [1, @ , whereas in 
SH-based models, hydrolysis occurs stochastically, any- 
where in the MT 0, H, El, [ll. Since IH models do not 
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allow for any GTP-bound tubulin within the GDP do- 
main, rescue depends solely on the probability of a free 
GTP tubulin attaching to a depolymersing MT with a 
GDP-tubulin at the tip. 

In SH models, GTP-monomers are located throughout 
the MT, with a concentration that decays exponentially 
with distance from the growing endp^. These models 
allow for a different rescue mechanism, i.e as the MT dc- 
polymerises by detachment of GDP-tubulins, GTP tubu- 
lins in the interior of the MT, or GTP 'remnants' get 
exposed and the MT starts polymerising again. Support 
for presence of GTP tubulins inside the MT has been 
provided by recent experiments Q. 

We study a particular SH model[l, [l^. In contrast to 
a SH model that takes into account all thirteen protofila- 
ments of a MT0, our model depicts the MT as a 1 — d se- 
quence with rates that prescribe polymerization, dcpoly- 
meriazation and hydrolysis. The appeal of this model is 
that there are only a few parameters, and the parameters 
can be controlled independently to understand the role 
that different processes play in the DI of MTs. 

The model represents the MT by a linear sequence 
of two species of monomers, which correspond to GTP- 
bound tubulin (denoted by -I- in rest of the paper) and 
GDP-bound tubulin (denoted by — ). We assume that the 
MT undergoes attachment and detachment only at one 
end, which we call the growing end (sometimes called the 
+ end in the literature). A MT evolves via the following 
rules: 

1. Attachment: If the growing tip is a + monomer, it 
grows with rate A by addition of a -I- subunit. 

2. Detachment: A — monomer at the growing end 
detaches with rate /i, causing its shrinkage. 

3. Hydrolysis: With rate 1 any -I- monomer in the MT 
can undergo hydrolysis to yield a — monomer. 

4. Non-remnant Rescue: + subunits could attach to 
a growing end with a — monomer at the tip with 
rate pX {p < 1). 

A previous study of the model [l^l for p > 0, demon- 
strated a transition from a phase of bounded to un- 
bounded growth of the MTs. The present study will focus 
on low p and fluctuations in the bounded growth region 
of the phase diagram. 
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Recent experiments monitored the distribution of 
lengths of growing and shortening excursions in in vitro 
systems of MTs. These experiments were able to re- 
solve fluctuations at the monomer level, and the distri- 
butions were found to be exponential. In simulations of 
our model we find that for small values of p the growth 
excursions arc independent of /i and can be fitted well 
by exp(— i/A)/A (i is the length of the excursion). The 
form of the distribution of shortening excursions changes 
non-trivially with jj. (Fig. [1]). 

We can fit the experimental data on the distribution 
of excursions [l^l with our model by taking A = 12; /x w 
0.1; <p< 0.001 (Fig. [J). The fitting provides tight 
bounds on the parameters, if we demand good fits to both 
positive and negative excursions simultaneously. For the 
low values of p obtained from the fits, rescue due to rem- 
nants dominates DI. Rare non-remnant rescues do in- 
fluence the steady state but have very little effect on the 
statistics of fluctuations. Hence, in order to highlight the 
effect of remnants on DI we analyze in detail the p = 
model, which we call the GTP remnants model since the 
only mode of rescue is via remnants. The presence of 
remnants and the correlation of their location with rescue 
events is consistent with recent experiments. In particu- 
lar, Perez et al. 0] observe GTP-bound tubulin within 
MTs and find that the location of these remnants corre- 
late to locations at which rescue events occurred during 
MT growth. 

In the remnants model, if the number of GTP 
monomers fluctuates to 0, rescue is no longer possible 
and the MT eventually shrinks completely. Any process 
that exposes remnants promotes rescues and makes the 
MT remain active for longer. In particular, increasing /i 
at fixed A leads to longer times (^at), and higher max- 
imum lengths (i(tjv)) of MTs before they disassemble 
due to loss of all GTP. 
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FIG. 1: Top: Comparison of excursions predicted by the GTP 
remnants model (x) with A = 12; ^ = 0.1, and experimental 
measurements (-I-) [l^. Bottom: The distribution of excur- 
sions with A = 12 for /i = 0.1, 0.5 and 1. While growth excur- 
sions have the distribution exp(— a::/A)/A, the distribution of 
shortening excursions changes from being exponential to non 
exponential and broader distributions with increasing /i. 



The time Im is a an interesting characteristic of DI in 
the GTP remnants model, since beyond Im the MT can 
only depolymerize and shrinks to zero in time of order of 
L[tiq) I fi. Fig[5]shows < tN > as a, function of In fj. for var- 
ious values of A. We see that < > increases monotoni- 
cally with increasing A and /i. For a given A as we increase 
fi, < t]^ > eventually saturates at a value of the order of 
exp(A)/A. In Fig. [21 we also show numerically obtained 
values of the growth velocity(w =< L{tN) > / < ti\r >) 
at time i jv • The velocity increases with increasing A and 
/i, and for a given A, it saturates for large /x. As illus- 
trated in Fig. 2, w = (A — + /.j) leads to good 
scaling collapse of the data. 

The distributions of and L{tisi) for various A and fi 
(Fig. [21) broaden with increasing ^, reaching asymptotic 
forms for /i >> A. The increase of fluctuations, indi- 
cated by these broadening distributions, is a consequence 
of remnants. Experimental measurements of these distri- 
butions should be able to distinguish between SH and IH 
mechanisms and provide indirect evidence for the pres- 
ence of remnants. Keeping this in mind, we provide some 
exact results for distribution functions in the limit of (a) 
/i = 0, and (b) A — > oo and /i — > oo. Obtaining approx- 
imate analytical expressions for distributions and aver- 
ages allows experimental verifications of the predictions 
of the model at different values of the system parameters. 

The probability distribution of total length of the 
MT(L(i)) at time t follows the following equation: 

= A(l - n,mP{L -l,t)- PiL, t)]+ 
^ino{t)[P{L + ^,t)-P{L,t)] 

where uq (t) is the probability of having a GDP at the tip 
of MT and is given by: 



dnojt) 
dt 



1 - noit) - nna{t)P{+- >,t) 



(2) 



here P(-| — >,t) is the probability that, given the tip of 
the MT is a — , the second last tubulin dimer is a -|-. Let 
us first consider fi = for arbitrary value of A. In this 
limit, no{t) = 1 — cxp(— t), and substituting in Eq. [21 
leads to: 



P{L,t) = e 



(-A(l- 



')) 



(A(l 



(3) 



The average length at any time < L{t) >= A(l — 
exp(— i)). Similarly, the average number of GTP-boimd 
tubulins at time t \s < T{t) >= Xtcxp{—t). Hence, the 
average time at which the amount of GTP goes to zero 
scales roughly as ln(A). This is different from the expo- 
nential dependence of on A for large fi indicated by 
the scaling collapse in Fig. [51 As seen in Fig. [21 < tjv > 
increase monotonically from a value of order O(lnA) to 
a value of order exp(A) as we increase fx from to oo. 

Another limit where we can obtain exact results for 
< L{t) > and < T{t) > is A ^ oo, fj, —>■ oo. In this limit. 
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In p. In H 

FIG. 2: Left: The average time at which the number of GTP-bound subunits goes to zero (< tjv >) as a function of In fi. The 
inset shows a scaling collapse of the plots, with y-axis scaled by e^/A and the a;-axis displaced by A, for A = 8, 10, 12. Right: 
The average growth velocity {v =< L{tN) > / < tN >) for indicated values of A as a function of In^. In the inset, we have 
scaled v by ^(A — l)/(/i + 1) to show the scaling collapse for A = 5, 8, 10, 12 and /i > 1. The dataset which does not fall on the 
same curve correponds to A = 3. 



P(H — >,t) can be assumed to be equal to unity for all 
t. We found numerically that most catastrophes were of 
0(1) and hence this approximation provides a sensible 
description of the MT behaviour in the growing phase 
for any A and ji. With P(H — >, i) = 1, we get. 



\ 1 + 1^1 J V 1 + 

where is the modified Bessel function of first kind 

and Eq. pleads to: 



< L(t) >-- 



(1 + Ai) 



(5) 
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FIG. 3: (a)-(c) The distributions of times at which the 
amount of GTP in the MT goes to zero(t]v). The x-axis 
is scaled by In A (the average value of tjv for = 0). The 
distribution shifts to the right with increasing until it sat- 
urates, (d) The distribution of maximum lengths for A = 5. 
The distributions of maximum lengths and tjv have similar 
dependencies on /i. 

The result that the MT is always in the growing 
phase for A > 1 is a consequence of the assumption 



P(H — >,i) = 1, which breaks down at some finite time 
for arbitrary values of A and /i. From Fig. O^Right), 
however, it is seen that the growing stage is very well 
approximated by Eq. [5] as the velocity obtained from 
simulations matches well with the expression of average 
velocity in Eq. [S] In fact, for a given t, the distribution 
of lengths (Eq. [4]) broadens with in a manner similar 
to that seen in the simulations (Fig. (3)). The variance of 
distribution of P{L,t) from Eq. Incomes out to be j^, 
which depends only on fi and t. 

Similarly, one can solve for the distribution and mean 
of T(i): 



< T{t) >= 



Xf-L Xexp{—{1 + fi)t) A(^ — l)exp{—t) 



1 



m(i + m) 



This equation also matches the simulation results, as we 
found that the average amount of GTP in the MT during 
the growing phase fluctuates around jj^^- 

A number of other predictions can be made in this 
limit of large A and fi, which applies to the growing 
phase, t << (<Ar). For example, the average number 
of catastrophes at time t, to leading order is given by 
< C{t) >= (^I'^^y} and the average number of GTP is- 
lands in the MT is predicted to be : . 

The predictions, presented above, reflect the sensitivity 
of dynamical properties to the depolymerization rate /z, 
and are a fingerprint of remnants. Experimental tests of 
these predictions can, therefore, provide insight into the 
nature of hydrolysis and rescue mechanisms in DI. 

The above analysis was restricted to p = in order 
to highlight the effect of remnants, and because the dis- 
tribution of growth and shrinkage excursions in in vitro 
experiments [lOl, could only be fitted with small values of 
p. As p is increased, the dynamics changes from rescues 
primarily due to remnants at small p to non-remnant res- 
cues that are also present in IH models at p ~ 1. Fig. [3] 
shows the time trace for MT length for a representative 
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FIG. 4: MT length L{t) (red lines) as a function of time for 
A = 5,At = 4 and p = (Left); p = 0.01{Right). The blue 
lines show the total amount of GTP as a function of time for 
the same runs. Arrows in the second plot indicate the rescues 
due to p > 0. 



run for p = and p = 0.01. Analysis of these trajecto- 
ries shows that a small, non-zero value of p introduces 
rare rescue events (indicated by arrows in the figure). 
These events change the overall length of MT, but the 
statistics of growth, rescue and catastrophes remain sim- 
ilar to p = 0. Measurements of these statistics should, 
therefore, provide tests of the remnants model of DI and 
provide clear distinction between IH models that do not 
have remnants and SH models which do. 

Our model can be easily extended to accomodate more 
detailed features of MTs while keeping the basic mech- 
anism of remnant-dominated rescues. For example, spa- 
tially varying hydrolysis rates due to the structure of 
MTS0 can be modeled by quenching some GTP-bound 
sites in our 1-d model. Similarly, the effect of motors that 
mechanically depolymerize MTs without dependence on 
GTP-states can be modeled by assuming a depoly- 
mcrizing rate p for both GTP and GDP-bound tubu- 
lins. Preliminary studies of the model with quenched 
disorder indicate that, although the time for which MT 
grows changes and there is a transition to unbounded 
growth function of percent of quenched sites, the 
distribution of excursions and velocity of growth remains 



unchanged, and therefore is a robust feature of the rem- 
nant model. Studies of the model mimicking motors also 
indicate that the basic features of the remnant model 
remain unchanged as long as A > /i -f 1 . 

To summarize, we have studied the role of rescue due to 
remnants in MT dynamics, and shown that remnants give 
rise to features of DI that are very different from IH mod- 
els. These features are robust, and with recent progress 
in experimental techniques d, , should provide tools 
for resolving the mechanism of hydrolysis inside a MT. 
Some particularly notable features of SH models, where 
remnants play an important role, that distinguish them 
from IH models are: 1) The distribution of MT lengths is 
not exponential but instead depends on A and fi, broad- 
ening as A and p, are increased. Measurements of rates 
of growth and shortening of indvidual MTs, both in-vivo 
and in-vitro, indeed exhibit variability that cannot be ex- 
plained by exponential distributions [15[. 2) The distri- 
bution of times for which a MT grows before a complete 
collapse also has a broad, non-exponential distribution in 
SH models. 3) The velocity of growth, besides depend- 
ing on the free tubulin conccntration( through A), also 
depends on depolymerisation and hydrolysis rates. Inter- 
estingly, at the same tubulin concentration, MTs exhibit 
much higher growth rates in-vivo in comparision to in- 
vitro [TtI . In conclusion, our study of a minimal model of 
MT DI, leads to predictions that provide sensitive tests 
of the mechanism of hydrolysis, and of the presence of 
GTP remnants in MTs. Experimental tests of these pre- 
dictions should lead to a better understanding of how 
DI is controlled through various agents both in-vivo and 
in-vitro. 
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